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Q\ ■ Abstract 

\ The spontaneous vacuum magnetization at finite temperature is inves- 

tigated in SU{2) gluodynamics within a consistent effective potential ap- 
^h. proach including the one-loop and the correlation correction contributions. 

<-j \ To evaluate the latter ones the high temperature limits of the polarization 

operators of charged and neutral gluon fields in a covariantly constant 
•H ■ magnetic field and at high temperature are calculated. The radiation mass 

squared of charged gluons is found to be positive. It is shown that the fer- 
romagnetic vacuum state having a field strength of order (gH) 1 / 2 ~ g^ 3 T 
is spontaneously generated at high temperature. The vacuum stability and 
some applications of the results obtained are discussed. 



1 Introduction 

The problem on generation of magnetic fields in nonabelian gauge theories at 
finite temperature is of great importance for particle physics and cosmology. Its 
positive solution, in particular, will give a theoretical basis for investigations of 
the QCD vacuum at high temperature and the primordial magnetic fields in the 
early universe |L|, [Q. In literature different mechanisms of producing the fields 
are discussed (see recent papers @,(§] and references therein). The aim of the 
present paper is to investigate in more detail one of them - the spontaneous 
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magnetization of the vacuum of non-abelian gauge fields at finite temperature. 
This problem was studied recently in Refs.|jl|, ||, || where the creation of the 
vacuum magnetic field has been derived. However, a trusty conclusion about 
the possibility of this phenomenon as well as estimations of the value of the 
condensed field have not been obtained, yet. In fact, in papers |jlj], the effective 
potential (EP) was calculated in one-loop order. But, as it is known, at finite 
temperature, T, the correlation corrections to the EP are of great importance. 
They may in an essential way influence the properties of a system. In Ref. H 
a number of the corrections has been taken into account. These authors have 
determined that the possibility of the vacuum magnetization depends on the 
sign of the next-to-leading term of the polarization operator (PO) of charged 
gluons, which has to be calculated in the field H and at T ^ 0. However, these 
calculations had not been done there. Besides, in Refs.[|l|,[|] the temperature 
mass squared of the unstable mode present in the spectrum of charged gluons 
was erroneously identified with the Debye one having the order ~ g 2 T 2 . However, 
the latter mass is generated for longitudinal modes of gauge fields. The unstable 
mode e 2 {n = 0, cr = 1) = p\ — gH, where g is the gauge coupling constant, 
n is the Landau level number, a is the spin projection on the field direction, 
is the transversal one produced by a spin interaction with the magnetic field. 
Therefore, the vacuum stabilization condition at T ^ has not been studied in 
detail. Moreover, a part of diagrams describing correlation corrections was missed 
there. In papers the scale of the magnetic field at T ^ was incorrectly 

estimated and as a consequence the authors came to negative conclusion about 
the vacuum magnetization. 

To summarize the present time situation we recall that in one-loop approx- 
imation the Savvidy level of order (gH))J 2 ~ g 2 T is generated. But the way of 
its stabilization and its existence with the correlation corrections been taken into 
account remain to be investigated. 

It will be important for what follows to remind recent results on observa- 
tion of the gluon magnetic mass in lattice simulations which was found to be 
of order m mag ~ g 2 T (as it has been expected from nonperturbative calculation 
in quantum field theory p(1,P0|). The mass screens magnetic fields at distances 
I > lm ~ {g % T)~ l but inside the space region I < l m they may exist. Since the 
typical values of particle masses at high temperature are M ~ gT, the magnetic 
fields of order (gH) 1 / 2 ~ gT or at least ~ g 2 T are to be of interest. They are 
able to affect all the processes at high temperatures. 

In the present paper we investigate the vacuum magnetization at finite tem- 
perature within SU(2) gluodynamics. Considering the Abelian covariantly con- 
stant chromomagnetic field H a = 5 a3 H = const and finite temperature as a 
background we calculate the EP containing the one-loop and the ring diagram 
contributions of both neutral and charged gluon fields. To find the latter ones the 
high temperature limits of the PO of gluon fields at the magnetic background are 
also calculated. It will be shown that in the adopted approximation the Savvidy 
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level with the field strength (gH)^ 2 ~ g 4 ^ 3 T is produced. Thus, we come to the 
conclusion that the correlation corrections make the field strength stronger as 
compare to the value derived in one-loop approximation This, in partic- 

ular, means that although the field is screened for I > (g 2 T) _1 , the spectrum of 
charged particles, being formed at Larmor's radius r ~ (gH)^ 1 ^ 2 ~ (g 4 / 3 T) -1 , is 
located inside this domain r «l m for small g. Moreover, since at high tempera- 
tures the particle masses are of order ~ gT , one has to conclude that a constant 
background magnetic field is a good approximation to investigate processes at 
high temperature. 

The content is as follows. In Sect. 2 the one-loop EP (zero and finite tem- 
perature parts) is presented in the form convenient for numerical investigations. 
Here we also present a general expression describing the contribution of the neu- 
tral gluon ring diagrams. In Sect. 3 we calculate the Debye mass of the neutral 
gauge field at H, T ^ 0. In Sect. 4 the same is done for the charged gluon part. 
Sect. 5 contains calculation of the two-loop diagram taking into account the self- 
interaction of charged gluon in the vacuum. The investigation of the vacuum 
magnetization and its stability as well as the discussion of the results obtained 
are presented in Sect. 6. In the APPENDIX we adduce the necessary high 
temperature asymptotic formulae. 



2 Basic Formulae and General Considerations 

Let us consider the pure Yang-Mills action 

S = -\ J d 4 xF« u Fr, (1) 

where 

f% = d »K ~ d « A l + 9f ahc A\K (2) 

is the field strength. In the SU(2) gauge theory the Savvidy vacuum state is 
characterized by the uniform background classical colour magnetic field 

Et = 5 a 3 5 l3 H, (3) 

where if is a constant, a = 1, 2, 3 is the isotopic index. 

To introduce this background field one has to decompose the gauge field po- 
tential as 

A a = A a ext + 3. (4) 

where ext = S^^Hxi and B^x) is the quantum field. In the EP calculations 
the background gauge fixing condition 

d,Bl^gr hc A h ^Bl = ^ (5) 

will be used. 
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The thermodynamic potential of the model is 



n = -h og z, (6) 

and 

Z = Trexp(-f3H) (7) 

is the partition function, 7i is the Hamiltonian of the system, /3 — 1/T is inverse 
temperature and the trace is calculated over all physical states. 

For what follows it will be convenient to introduce the "charged basis" of 
fields 

Wf = ±(Al±iAl), A, = A% (8) 

In this basis the problem of the vacuum magnetization is reduced to calculation 
and investigation of the vacuum polarization of fields W^, W~ in the external 
field ,l' w . 

To obtain the EP one has to rewrite eq.(|7]) as a sum over quantum states 
calculated near the nontrivial classical solution A ext . This standard procedure is 
described in many papers and textbooks (see, for example, Refs. [ . [ 1 1 1 . [|13H . 
H) and the result can be written in the form: 

V = -~ £ TrlogG? ) + V&(H, T) + • • • , (9) 
z i 

where i marks the type of field (W ± ,A and ghosts), is the correspond- 
ing propagator in the external field A ext . In the present paper to incorporate 
temperature the imaginary time formalism will be used. In the specified above 
background field, the trace means summation over the discrete Matsubara fre- 
quencies, summation and integration over eigen values of the quantum fields. 
The first term in Eq.^) corresponds to the one-loop EP whereas the other ones 
present the contributions of two-, three-, etc. loop corrections. 

Among these terms there are ones responsible for dominant contributions of 
large distances at high temperature - so-called daisy or ring diagrams (see, for 
example Ref. |fI2| ). This part of the EP, V ring (H,T), is important in the case 
when massless states appear in a system. The ring diagrams of either neutral 
and charged gluons have to be calculated when the vacuum magnetization at 
finite temperature is investigated. Really, one first must assume that the field is 
nonzero, calculate EP V(H, T) and after that check whether the minimum of it 
is located at nonzero H. On the other hand, if one investigates problems in the 
applied external field, the charged fields become massive with the mass ~ (gH) 1 ^ 2 
and have to be omitted. To find V r i ng (H,T) the one-loop polarization operators 
of charged, II(if, T), and neutral, U°(H,T), gluons in the external field and at 
finite temperature have to be calculated in the limit of zero momenta. Then 
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V r i n g(H,T) is given by the series depicted in figures 1,2. Here, the dashed lines 
describe the neutral gluons and the wavy lines represent the charged ones, blobs 
stand for the one-loop polarization operators. The diagrams with one blob show 
the two-loop terms of the EP, with two blobs - three-loop ones, etc. 




Figure 1: The neutral gluon ring diagrams giving contribution to the effective 
potential. 





Figure 2: The charged gluon ring diagrams giving contribution to the effective 
potential. 

The standard procedure to incorporate these contribution is to substitute df^ 
in Eq. (|) by the inverse full propagator Gj l = [Gf 5 ] -1 + U(H,T) p], @. In 



the case under consideration the expression for the EP can be presented as follows 
= £ J^^tjog[P^f + el am +U(n,T,H,a))] (10) 



oo 
oo 



n=0,o-=±l 



W i^l (2-) 

This generalized EP is written as the sum of energies of the charged gluon field 
modes in the external magnetic field ( uji — ^ - discrete imaginary energies) 

el a =pl + (2n + l-2a)gH. (11) 

It includes the temperature masses U(H,n,cr,T) of the charged modes, which 
are also dependent on H, the level number n = 0, 1, ... and the spin projection 
a = ±1, as well as the one of the neutral gluons, U°(H,T). Then, denoting the 
sum of + e 2 describing the tree level spectrum as Dq 1 ^, H,T), let us rewrite 
the first term of Eq.flTUD as follows 
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°°-oo n,a 



+ m £ /°^{io g [i + ( w f+^-^- 1 n(^T)] 



gH_ f dp 3 
£ log[l + L> (e^F,T)n(F,r,n )( 7)]}. (12) 



Here, the first term gives the one-loop contribution of the charged gluons, 
calculated and investigated in detail in Refs. U , H , the second one is the sum 
of the ring diagrams of the unstable mode e^ =0 cr=+1 = pl — gH (as it can be easily 
checked by expanding the logarithm in a series) and the last term describes the 
sum of ring diagrams of the stable modes. 

In general, the incorporation of the polarization operators into EP may lead 
to incorrect account for the combinatoric factors of the two-loop diagrams. In our 
case, to be sure in the correctness of the obtaining results the following procedure 
has been applied. We subtracted the term 

v s = £ / ^ £ d °(p*> H > T ) n (^ n > T > ^ ( 13 ) 

l — OO — 



n = Q, 

CT = ±1 



from Eq. (p!Q|) that separates the contributions of the two- loop diagrams of charged 
gluons from correlation corrections (the first diagram in Fig. 2). Then, it is easy 
to check that the two-loop diagrams containing the one of loops having the neu- 
tral gluon lines (the first diagram in Fig.l) can be accounted for within the 
second term of the EP ([TD|) . After that the only two-loop diagram describing the 
self-interaction of charged gluon fields in the vacuum has been straightforwardly 
computed. 

In one-loop order the neutral gluon contribution is a trivial if-independent 
constant which can be omitted. However, these fields are long-range states and 
they do give if-dependent EP through the correlation corrections depending on 
the temperature and field. Below, only the longitudinal neutral modes are in- 
cluded because their Debye's masses U°(H,T) are nonzero. The corresponding 
EP is easily calculated and is given by the expression which can be recognized 



from H = case HI, [113 



Vnn 9 = ^U°(H,T)T 2 -^-[U°(H,T)f 2 

(Jl°(H,T)) 2 n . 4ttT , 3 
+ g 27r2 (nop) + 4 -7], ( 14 ) 
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and H°(H, T) = n(j (A; = 0, H, T) is the zero-zero component of the neutral gluon 
field polarization operator calculated in the external field at finite temperature 
and taken at zero momentum, 7 is Euler's constant. [] The first term in Eq. (|T4|) 
has the order ~ g 2 in the coupling constant, the second term is of order ~ g 3 and 
the last one - ~ g 4 . Restricting ourselves by the order ~ g 3 it will be omitted in 
what follows. As usually, for T1(H, T) the high temperature limits of the functions 
have to be substituted into Eqs. (|14|), dl2|)[^|, As we noted above and one 
is able to check, this part of the EP correctly accounts for the contributions of 
the two-loop diagrams with the neutral gluon lines. 

The detailed calculations of the one- loop EP (H, T) have been carried out 
in Refs. ||, @. The EP can be written in the form: 

v m = J_ 7% e -^ g e^(iiVM[ gHm ff, s) - -l (15) 

8tt 2 7 s 2 L sm(gHs) s 

where an auxiliary mass parameter /i 2 — ie, e — ► 0, is introduced which regulates 
an infrared region and plays the role of the normalization point in the field. It is 
useful for analytic continuations from weak fields gH < fi 2 to the field strengths 
gH > n 2 when an imaginary part of the EP is calculated. Then it can be set to 
zero. The term with / = gives the well studied zero temperature part, while 
other terms describe the statistical part. The proper time (s- representation) is 
used. By integrating over s we present Eq.(|l5D in the form convenient for numeric 
investigations M: 

V {1 \H,T) = V W (H) + V T W (H,T), 

w (H) = ^ + w!|li logW , v ^], 

71 P 1=1 p =0 1 1 1 

- ^i(^) 1/2 )]-i^i^Ey^(^(^) 1/2 ). (16) 

where Ki(x), Y±(x) and J\{x) are Bessel's functions. Remind that \i is a sub- 
traction point in the field H . The imaginary part of the EP is signaling the 
vacuum instability. To have a more transparent expression let us calculate the 
high temperature limit T >> (gH) 1 / 2 >> /x of the V^(H,T). By means of the 
Mellin transformation technique we obtain 0,0: 

„ m . rrm , H 2 llg 2 TT2 , T 2 l(gHf' 2 T 
2 48r ^6 71 

- i ( -^f^ + 0(g 2 H 2 ). (17) 

Z7T 



1 This type of the EP was used in Ref. Q to account for the correlation corrections of charged 
gluons. However, the corresponding correct expression is quite different (see Eq.(p9|)). 
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Note the important cancellation of the if-dependent logarithms entering the vac- 
uum and the statistical parts. 

Now, let us consider the second term in Eq.(^). An elementary integration 
gives 

VunstaUe = ^^(H, T, U = 0, O = +1) - gH] 1 ' 2 + JiMF-E. . (1 8 ) 

From Eqs. ([PTD and (18) it is seen that imaginary terms are cancelled out in the 
total. The final EP is real if the relation U unsta bie{H,T) > gH holds. To check is 
it the case or not one has to calculate both the field spontaneously generated in 
the vacuum and the radiation masses of charged gluons. The former value can be 
immediately calculated by differentiating Eq.(|I7|) with respect to H and setting 
the result to zero |J,|§: 



{gHfJ 2 = f^T. (19) 

To answer the question whether this vacuum is stable or not one has to com- 
pute U(H,T,n,cr), which is the average value of the gluon PO taken in the 
tree level state of the charged gluon, and consider the effective mass squared 
M 2 (H C , T) = H(H C , T, n = 0, a = +1) —gH c . If it is positive, one has to conclude 
that radiation corrections stabilize the vacuum (in one-loop order), otherwise it 
remains unstable. 



3 Debye Mass of Neutral Gluons 

First, let us calculate the if -dependent Debye mass of the neutral gluons. The 
following procedure will be applied. We calculate the one-loop EP due to vacuum 
polarization of charged gluons in the external field H and some chemical potential, 
p, which plays the role of an auxiliary parameter. Since the EP is the generating 
functional of the one-particle irreducible Green's functions of the field p, one has 
by differentiating it twice with respect to p and setting p = to obtain the mass 
squared m 2 D neutral = U 00 {H,T,p = 0). 

The corresponding temperature dependent part of EP is 

y (1) = "f^E J ^ exp HW4*)[— E ^y + 2sinh(^)]cosh(^p). (20) 

All the notations are obvious. Then, after the differentiations we get 



m D 



d 2 V^ gH d » [da l 2 (3 2 eH ^ 1 

~^~ip=o = S / T exp( -^^ } [ ^Hs) + 2 smh(s)] " 

(21) 
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Expanding sinh 1 s in series over Bernoulli's polynomials, 



1 e~ s ^ B k 



Elft- 2 ^' (22) 



sinhs s k\ 

and carrying out integration over s in accordance with the standard formula 

[ dss* 1 - 1 exp(-as - -) = 2(-) n/2 K n (2\^b), (23) 
J s a 

o 

a, b > 0, we obtain in the high temperature limit (gH) 1 / 2 /T « 1, 

i» = f /? 2 | D8 *ggy + 4AV( 9 ff)'^) + 0(ffl]- (24) 



Hence, summing up series by means of Mellin's transformation (see the Appendix) 
and differentiating with respect to /? 2 , one finds the Debye mass of neutral gluons, 

ml = \g 2 T 2 - - -L^U) + 0{{gHf/T 2 ). (25) 

3 7r 47r z 

Here, the first term is the well known temperature mass squared of gluon and 
other ones give field-dependent contributions. As it is seen, they are negative. 
This is important for what follows. It is also interesting that spin does not 
contribute to the Debye mass in the leading order. Substituting expression (f2~5|) 
into equation flI3|), we obtain the correlation corrections due to neutral gluons. 

4 Ring Diagrams for Charged Gluons 

Now, we are going to calculate U unsta bie(H, T) and H(H, T, n, a) which cannot be 
found from any effective potential and require an explicit calculation of the mass 
operator of the charged gluon field. This is a separate and sufficiently complicate 
problem which is discussed in detail in the other publication [15|]. Here, we only 



adduce the necessary for our present purpose results - the high temperature limits 

Of ^unstable (H,T): 

Tl unstable (H,T) =< n = 0,(7 = 1 | Relif u ar9ed \ n = 0,a = 1 >= Ibm^gHf^T, 

(26) 

and of excited states U(n ^ 0, a ^ +1), 

a 2 

ReU(p 4 = 0,n,p 3 = 0,H,T,a = +1) = f-(##) 1/2 T(15.62 + An), (27) 

47T 

a 2 

Reliipi = 0,n,p 3 = 0,H,T,a = -1) = — (gH) 1/2 T(UM + An), 

An 
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where the average values of the PO in the states of the spectrum ([11]) are cal- 
culated. These formulae correspond to the limit gH/T 2 « 1. The operator 
contains also an imaginary part which describes the decay of the states due to 
the transitions to lower energy levels. But for the problem under consideration 
only the real part is needed, since it describes the radiation mass. 

Let us note the most important features of the expressions (|26| ) ,(|27| ) . It is 
seen, at if = no screening magnetic mass is produced in one-loop order, as 
it should ||. Second, the masses squared of the modes are positive and act to 
stabilize the spectrum of charged gluons at high temperatures. Thus, one has to 
conclude that in a nonzero chromomagnetic field the charged transversal gluons 
become massive at finite temperature. (] 



The Debye mass of charged gluons is found to be TT5 



Tl o(h = 0, h = 0, H, T) = \g 2 T 2 + ^-( g E) x l 2 T{§ + An), (28) 

where again only the real part is presented. As it is seen, the next-to-leading 
term is the positive growing function of n. It is interesting to note that the 
Il(p4 = 0, P3 = 0, n, a, H, T) is more informative as compared the function 11(^4 = 
0,p = 0,cr, T) used at zero field. This is because it exactly accounts for the 
transversal momenta described by n variable: p\ + p\ — (2ti + l)gH. 

Now, substituting the the expressions (P6),(P7|), fl2"5p into Eq.flT2"D and inte- 



grating over momentum and calculating the sums in n, we obtain the ring diagram 
contribution of charged gluon fields. The result can be expressed in terms of the 
generalized ^-function and looks as follows, 

V r f n9 = ^{^/^[C(-\, a + ) +((—, a-) +2a~ a D )} 

+ (U(H,T,n = 0,a = +l)) 1 / 2 }, (29) 

where the first term in the first squared brackets corresponds to the spin pro- 
jection a = +1, the second term - a — — 1 and the last one describes the part 
due to longitudinal charged gluons. The terms in the second squared brack- 
ets give the independent of H(H, T) part of eq. (|T2"D . The last term in the 
curly brackets is due to the radiation mass of the unstable mode. This ex- 
pression is real for sufficiently high temperatures. The notations are introduced: 
9 H D =9 H+ i{ 9 H)WT , a_ = I + a + = \ 4- and 

°o - (jg 2 r 2 + g(gg) 1/2 r + g_H )/2 g H D . 

2 This conclusion is in obvious contradiction with the result of Refs. , Q where no radiation 
mass for the state with n — 0, a = +1 has been derived. These authors have included the 
dependence of H(n = 0, a = +1, H, T) on H through the eigen states of the tree-level spectrum, 
only. As the operator they used the zero field expression in the high temperature limit. But 
for transversal modes at P3 = this is zero |9j . 
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Substituting the expression ( p5|) into eq.(|T4|) and gathering all other contri- 
butions ([16|) , ([2~9"D, we obtain the consistent expression for the EP. To correctly 
account for the two-loop diagram contributions we have to subtract from it the 
terms (0), which will be calculated below. 



5 Two-Loop Contribution of Charged Gluons 

In this section, to complete the calculation of the EP, the contribution of the two- 
loop vacuum diagram taking into account the self-interaction of charged gluons 
is computed. As before, we use the Feynmann gauge and the Furry picture for 
Green's functions in the coordinate representation. If one denotes the propagator 
of the charged gluons in the external field H as iG^ u (x,y), the contribution of 
the diagram to the EP can be written as 

2 Q 2 

v ch = - — [G IItl (x,x)G uu (x,x)+G lliy (x,x)G Uf ,(x,x)-2G flu (x J x)G flu (x,x)], (30) 
where the Green function is [II 



I y I 

G tlu {x,y) = -—exjp{-ig J dz[A(z) + -F ■ (z — y)]} (31) 

X 

OO 

r ds r 1, sinhgFs i, . onF* 

J ^2 ex Pi~ 2 g ~F~s 4^ ~ V ' ' 9 Fcoth 9 Fs -\x-y) }e^ 9 . 

o a 

Here, A^(x) and are the potential and the field strength tensor of the ex- 
ternal field. The proper time ( s-representation) is introduced. To incorporate 
temperature in Eq. ( ftip the method of Ref. |R| will be used. After rotation to 
the Euclidean space the Matsubara-Green function is expressed in terms of the 
Green function at zero temperature by the formula 

+oo 

G(x, y-T)= (-l) n+W)A G(x - [x]/3u, y - n/3u), (32) 

n=~ oo 

where the corresponding Green function at T = enters the right-hand side, 
(3 — k,u — (0, 0, 0, l),[x] denotes the integer part of x±f j3 and we use the 
parameter A = 1 for fermions, A = for boson and ghost fields. 

Substituting Eq. (|32|) , where G in the right-hand side is given by Eq. (pT|) , into 
expression (|30|), we obtain the contribution V}£\H,T). By calculating the trace 
and setting y = x in G flu (x,y,T) the temperature dependent part of it can be 
written as follows, 

2 oo 

V%\H, T) = ^ £ [II (n, H, (3) + I 2 (n, H, (3) - 2I 3 (n, H, f3)), (33) 



11 



where the notation are introduced: 



h(n,H,(3) = 2iJ 



h - 2U 



+ 2 



o 

oo 



ds gH n 2 ,a 2 , 

e 4s (1 + cosh 2g Hs) 

s smhgHs 



ds gH i2 

6 4s 

s sinhgHs- 



dS n 2 /3 2 



-e 4s gHcothgHs 



dS n 2 (3 2 



-e 4s gH 



(34) 



The zero temperature part is given by the term with n 
discrete frequencies 
malization 



fsee Ref.[16 



in the sum over the 
It contains divergences and require a renor- 



The statistical part is finite. 
To find the high temperature limit of these expressions we first carry out 
the integration over the parameter s, as in Sect. 3. That can easily be done by 
expanding sinh -1 (gHs) and coth(gHs) in series over the Bernoulli polynomials 
and integrating over s in accordance with formula (|23|). The calculation of sums 
in n can be carried out by means of the Mellin transformation technique (see, for 
example, Ref.|J). The necessary formulae are adduced in the Appendix. In this 
way we obtain for the leading terms of the limit of interest: 



V, 



S\H,T) lT _ = -l-(gH) 1 /^ + O(gHT') 



The negative sign is important for what follows. 
The high temperature limit of the term (O) is 



(35) 



V s (H,T) lT _ = ^=-C(^,-^)(^) 1/2 T 3 + 0(^T 2 ). 



(36) 



It was also calculated by the Mellin transformation method. This term must be 
subtracted from the total EP and therefore gives the negative contribution to 
the leading terms of the asymptotic expansion because £(f > — |) = 0.8093. It is 
worth to mention that these are the longitudinal modes that determine the high 
temperature behaviour of V^(H,T). Having obtained the two-loop corrections 
to the EP, one is able to investigate the spontaneous magnetization of the vacuum 
at high temperature. 



6 Discussion 

The derived EP is expressed in terms of the well known special functions. There- 
fore, it can easily be investigated numerically for any range of parameters en- 
tering. As usually, it is convenient to introduce the dimensionless variables: the 
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field = (gH) x ' 2 /T and the EP v((j),g) = V(H,T)/T 4 . The vacuum magnetiza- 
tion at high temperatures, T » (gH) 1 ^ 2 , 0^0, will be investigated within the 
following limiting form of the EP: 



11 4 , ,T 2 , 10 3 g 2 

■ total U n\>, „ = — 1 _ loef — ) — < 

2g 2 48 7T 2 *% 2 3 7T 24tt 

£ 3 1 # 2 , # 2 



■ 0,80930 + 0{g 3 ). (37) 



3^0+^ 4?r QV^tt 

The logarithmic term is signalling the asymptotic freedom of g 2 {T) at high tem- 
peratures |J. It includes explicitly the dependence on the scale parameter 
fi. Other terms present, respectively, the high temperature asymptotics of the 
one-loop EP, the neutral gluon and the charged gluon rings and the two-loop di- 
agram contributions. To obtain the term due to V^L the asymptotic expansion 
of the Zeta-function [^(J 



C(~2> a D)\a D ^oo = -\ a T + a D 2 - -^ a D /2 + °( a D /2 ) (38) 

has been used. Zeta- functions with a + ,a_ do not contribute in leading order. 
Since we are searching for the fields of order greater then g 2 , we can omit the 
term g 3 and obtain for the condensed field 

{9H)T = % 9 i/3 T. (39) 

Thus, we come to the conclusion that the ferromagnetic vacuum state does exist 
at high temperatures. The correlation corrections increase the field strength as 
compared to the value (gH) l J 2 ~ g 2 T generated in one-loop order. 

Now, let us compare the results of investigations in Refs. |1], [|J, with that 
of the present paper. In Ref . the one- loop EP was calculated. In equation ( |37|) 
this corresponds to the first and third terms which determine the vacuum mag- 
netization of order (gH) l J 2 ~ g 2 T.The authors of Ref. QXJ] missed the important 
cancellation of the logarithmically dependent on H terms presenting in both the 
vacuum and the statistical parts of the EP. Hence, their estimate of the magnetic 
field strength generated at the GUT scale is actually based on the zero temper- 
ature EP calculated by Savvidy fIT |. As one can see, there are no dependent 



on log(gH) terms in Eq. (|37|) . Further, in Ref. [0J the Debye mass ~ g 2 T 2 as 
the temperature dependent radiation one of the transversal charged gluons has 
been substituted, that is an incorrect assumption (although, in principle, it is not 
important from the point of view of the result of our investigation because the 
radiation mass of charged transversal modes has been found to stabilize the spec- 
trum at high temperature). The Debye mass for the transversal modes was also 
used in Ref. [[|. Therein the mentioned cancellation of the logarithms was estab- 
lished and in one-loop approximation the condensed field of order (gH)^ 2 ~ g 2 T 
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has been obtained. However, the temperature mass m 2 D was used twice. First, as 
a heuristic factor providing the vacuum stabilization at high temperature. Sec- 
ond, in the structure of the EP responsible for the ring diagram contribution of 
the unstable mode. The latter is the main point for the final result. Besides, the 
contribution of the neutral gluons was missed and the dependence on the Landau 
level number, n, has not been assumed at all. Therefore, the EP used in Ref.|| 
gave no possibility to calculate the correct value of the condensed field. 

The speculations on impossibility of the spontaneous vacuum magnetization 
in Refs. [0] , f| are based on the results of Ref . Q and the observations of the gluon 
magnetic mass at finite temperature of order ~ g 2 T. Since the field at finite tem- 
perature was identified with the zero temperature value (gif) 1,/2 (T = 0) << g 2 T 
whereas the spectrum of the charged particles is formed at distances / ~ (gH)~ l l 2 
much longer then the screening magnetic length the absence of the magnetization 
has been claimed. Q Our EP is free from these short comings. It includes all the 
relevant terms in leading order and the temperature masses with the correct de- 
pendences on H and T for all modes. With these improvements made we come to 
the conclusion that the ferromagnetic vacuum state is spontaneously generated at 
high temperature and the field strength is found to be of order (gH) l J 2 ~ g 4 ^ 3 T. 
For small g{T) this is stronger then the field determined in one-loop order. 

To better understand the situation with the vacuum stability, let us first con- 
sider the one-loop case. As the effective mass squared we find for the condensed 
field {gHfJ 2 = (g 2 /2n)T: M 2 (H,T) = U(H c ,T,n = 0, a = +1) - gH c > 0. 
Thus, the vacuum stabilization is observed. But if one substitutes instead this 
value the field strengths of order (gH)^ 2 ~ g 4 ^ 3 T described by equation (|3^) 
the effective mass squared becomes negative. The gluon one-loop radiation mass 
does not stabilize the true vacuum magnetic field. 

Nevertheless, the one-loop result makes hopeful the idea to have the stable 
vacuum due to radiation corrections to the charged gluon spectrum. Naturally, 
to investigate this possibility the gluon polarization operator with the correlation 
correction included should be calculated. This problem requires an additional 
investigation. Other interesting possibility is the formation at high temperatures 
of the gluon electrostatic potential, so-called A condensate (see survey pTf), 



which also acts as a stabilizing factor ||. To realize the latter scenario consistently 
the simultaneous spontaneous generation of both the Aq condensate and the 
magnetic field should be investigated. If again the homogeneous vacuum field 
will be found to be unstable with these improvements made, the inhomogeneous 
fields of the lattice type discussed in papers [§], |p2fl ,p3| may be created. Since 
the condensed magnetic field is strong at high temperatures, the lattice structures 
having the cells of order ~ l/(g 4 ^T) « l/{g 2 T) are located inside the region 

3 This estimate is used in Ref. jl8| devoted to investigation of cosmic magnetic fields in 
inflationary universe. However, in connection with the results presented above some of numeric 
estimations has to be corrected. 
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where the fields are not screened by the gluon magnetic mass. 

These results may found applications in problems of cosmology, in particular, 
in studying the primordial magnetic fields in the early universe. In a few recent 
years the scenario of the evolution of the universe in external hypermagnetic field 
[|J, |23| became popular. Our calculation of the vacuum magnetization unam- 
biguously determines the possibility of the presence of strong magnetic fields in 
the hot universe. Other application is the high temperature QCD. Here, one 
also should take into account the generation of the chromomagnetic field in the 
deconfinig phase. As a conclusion we stress again that Savvidy's mechanism of 
the vacuum magnetization does work at high temperature, although a number 
of questions concerning the vacuum stability has to be investigated in order to 
derive a final picture. 

The authors thank J. Bohacik for useful discussion. This work has been partially 
supported by DFG grant No. 436 UKR 17/24/98. 



Appendix 

To calculate the the sums in n in Sects. 3, 5 the Mellin transformation was 



used. It is described in detail in Refs. flSSLKJ. Here, we adduce the results for 



summations of series appeared in the considered cases. In high temperature limit 
gH/T 2 « 1 ( uj — > )one obtains: 

£ -K^uuj) = C -^~l + j(log - - 7 + 1/2) + 0(u, 2 ), (40) 

°° 1 47T 

£ K (nuj) = -( 7 - log -) + 0(u), (41) 

n=l 1 U 

where 7 is Euler's constant. 
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